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Our  research  for  1990  has  been  widely  reported,  as  given  by  the  listing  following,  of  16  Journal 
Articles,  3  ERL  Reports,  4  Talks,  and  12  Poster  Papers. 

Abstracts  are  attached  for  some  of  the  talks. 

Sent  along  with  this  Report  are  reprints  of  Journal  Articles  and  the  ERL  Reports. 

Our  prior  mode  was  to  publish  Quaterly  Progress  Reports;  these  then  became  Semi-Annual 
Reports,  which  ended  in  1988.  In  1989,  we  began  publishing  Annual  Progress  Reports.  While  QPR’s 
were  excellent  exercises  in  reporting,  they  required  an  immense  effort;  in  today’s  research  climate,  such 
effort  is  not  available. 

We  trust  that  our  reporting  is  still  useful!. 

C.K.  Birdsall 
Principal  Investigator 


Publications  for  1990 


Journal  Articles 


M.  Surendra,  D.B.  Graves,  and  IJ.  Morey,  "Electron  Heating  in  Low-  Pressure  RF  Glow  Discharges,” 
Appl.  Phys.  Lett.,  56,  pp.  1022-1024,  March  12,  1990. 

Lou  Ann  Schwager  and  C.K.  Birdsall,  "Collector  and  Source  Sheaths  of  a  Finite  Ion  Temperature 
Plasma,"  Phys.  Fluids  B2,  pp.  1057-1068,  May  1990. 

MJ.  Gerver,  S.E.  Parker,  and  K.  Theilhaber,  "Analytic  Solutions  and  Particle  Simulations  of 
Cross-Field  Plasma  Sheaths,"  Phys.  Fluids  B2  pp.  1069-1082,  May  1990. 

IJ.  Morey  and  C.K.  Birdsall,  "Traveling  Wave-Tube  Simulation:  the  IBC  Code,"  IEEE  Trans.  Plasma 
Sciences.,  pp.  482-489,  June  1990. 

J.K.  Lee  and  C.K.  Birdsall,  "Particle  Simulation  of  Drift  Cyclotron  Instability",  Japanese  Journal  of 
Applied  Physics,  Part  1,  29,  pp.  2126-2134,  October  1990. 

RJ.  Piocassini,  C.K.  Birdsall,  B.I.  Cohen,  "Particle  Simulation  of  Collisional  Transport  in  a  Diverted 
Scrape-off  Layer.  Part  I-The  High-Recycling  Regime, 'Wuc/.  Fusion.  30  pp.  2329-2347,  November 
1990 

RJ.  Procassini,  C.K.  Birdsall,  E.C.  Morse,  "A  Fully-Kinetic,  Self-Consistent  Particle  Simulation  Model 
of  the  Collisionless  Plasma-Sheath  Region,"  Phys.  Fluids  B2  pp.  3191-3205,  December  1990. 

S.E.  Parker,  and  C.K.  Birdsall,  "Numerical  Error  in  Electron  Orbits  with  Large  0)„  At ,"  accepted  by  J. 
Comp.  Physics. 

A.  Friedman,  S.E.  Parker,  S.L.  Ray,  and  C.K.  Birdsall,  "Multi-Scale  Particle-In-Cell  Plasma  Simula¬ 
tion,"  accepted  by  J.  Comp.  Physics. 

S.E.  Parker,  A.  Friedman,  SJL.  Ray,  and  C.K.  Birdsall,  "Bounded  Multi-Scale  Plasma  Simulation: 
Application  to  Sheath  Problems,"  accepted  by  J.  Comp.  Physics. 

N. F.  Otani,  J-S  Kim,  C.K.  Birdsall,  B.I.  Cohen,  W.  Nevins,  N.  Maron,  "Elimination  of  Velocity  Space 
Rings-and-Spokes  Instabilities  in  Magnetized  Electrostatic  Particle  Simulations  of  Plasmas,"  accepted 
by  J.  Comp.  Physics. 

S.E.  Parker,  RJ.  Procassini,  C.K.  Birdsall,  B.I.  Cohen,  "A  Suitable  Boundary  Condition  for  Bounded 
Plasma  Simulation  without  Sheath  Resolution,"  accepted  by  J.  Comp.  Physics. 

V.  Vahedi,  M.A.  Lieberman,  M.V.  Alves,  J.P.  Verboncoeur,  and  C.K.  Birdsall,  "A  One  Dimensional 
Collisional  Model  for  Plasma  Immersion  Ion  Implantation,"  accepted  by  J.  Appl.  Physics. 

M.V.  Alves,  M.A.  Lieberman,  V.  Vahedi,  and  C.K.  Birdsall,  "Sheath  Voltage  Ratio  for  Asymmetric  RF 
Discharges,"  accepted  by  J.  Appl.  Physics. 

J.P.  Verboncoeur,  M.V.  Alves,  and  V.  Vahedi,  "Simultaneous  Potential  and  Circuit  Solution  for 
Bounded  Plasma  Particle  Simulation  Codes,"  submitted  to  J.  Comp.  Physics  June  1990. 
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T.L.  Crystal,  P.C.  Gray,  W.S.  Lawson,  C.K.  Birdsall,  S.  Kuhn,  "Trapped  Electron  Effects  on 
Time-Independent  Negative-Bias  States  of  a  Collisionless  Single-Emitter  Plasma  Device:  Theory  and 
Simulation,"  to  appear  in  Physics  of  Fluids  B,  January  1991. 


ERL  Reports 


M.V.  Alves,  M.A.  Lieberman,  V.  Vahedi,  and  C.K.  Birdsall,  "Sheath  Voltage  Ratio  for  Asymmetric  RF 
Discharges,"  Memo.  No.  UCB/ERL  M9Q/56,  June  21, 1990. 

V.  Vahedi,  M.A.  Lieberman,  M.A.  Alves,  J.P.  Verboncoeur,  and  C.K.  Birdsall,  "A  One  Dimensional 
Collisional  Model  for  Plasma  Immersion  Ion  Implantation,"  Memo.  No.  UCB/ERL  M90/60,  July  9, 
1990. 

J.P.  Verboncoeur  and  V.  Vahedi,  "PDP1:  Plasma  Device  Planar  1  Dimensional  Bounded  Electrostatic 
Code.  Reference  Manual,"  August  1990. 


Conference  Proceedings,  Poster  Papers 


Sherwood  Fusion  Theory  Conference,  Williamsburg,  VA,  April  23-25,  1990: 

RJ.  Procassini,  and  C.K.  Birdsall,  "Particle  Simulations  of  Transport  in  a  Diverted  Tokamak 
Scrape-Off  Layer  The  High-Recycling  Regime." 

Live  demonstration  in  Plasma  Visualization  evening  session  7- 11pm  April  24  of  our  bounded 
plasma  PC  codes,  applications,  (by  Birdsall,  Vahedi) 

Microwave  Power  Tube  Conference,  Monterey,  CA,  May  7-9,  1990: 

IJ.  Morey,  and  C.K.  Birdsall,  "Traveling-Wave  Tube  Simulation;  the  IBC  Code." 

1990  IEEE  International  Conference  on  Plasma  Science,  Oakland,  CA,  May  21-23,  1990: 

R  J.  Procassini,  C.K.  Birdsall,  B.I.  Cohen,  "Particle  Simulations  of  Transport  in  a  High-Recycling 
Divertor  Scrape-Off  Layer." 

M.V.  Alves,  M.A.  Lieberman,  V.  Vahedi,  C.K.  Birdsall,  "Sheath  Voltage  Ratio  For  Asymmetric 
RF  Discharges." 

J.P.  Verboncoeur,  V.  Vahedi,  M.A.  Lieberman,  C.K.  Birdsall,  "Work  Done  And  Energy  Balance 
in  RF  Discharges." 

V.  Vahedi,  M.A.  Lieberman,  M.A.  Alves,  J.  P.  Verboncoeur,  C.K.  Birdsall,  "A  Collisional  Model 
For  Plasma  Immersion  Ion  Implantation." 

Workshop  on  High-Density  Plasma  Techniques  and  Processes  for  Integrated  Circuit  Fabrication,  Bur¬ 
lingame,  CA,  September  11-12, 1990. 

V.  Vahedi,  M.A.  Lieberman,  M.A.  Alves,  J.P.  Verboncoeur,  and  C.K.  Birdsall,  "A  One  Dimen¬ 
sional  Collisional  Model  For  Plasma  Immersion  Ion  Implantation." 

V.  Vahedi,  J.P.  Verboncoeur,  M.A.  Alves,  and  C.K.  Birdsall,  "Plasma  Processing  Via  Computer 
Simulation." 
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43rd  Annual  Gaseous  Electronics  Conference,  Urbana-Champaign,  IL,  October  16-19,  1990. 

C.K.  Birdsall,  "Particle -In-Cell  Combined  with  Monte  Carlo  Collisions  In  Living  Color."  (Invited 
talk) 

J.P.  Verboncoeur,  V.  Vahedi,  and  C.K.  Birdsall,  "Power  Deposition  in  Parallel  Plate  Discharges." 

APS/Division  of  Plasma  Physics,  Cincinnati,  OH,  November  12-16,  1990. 

S.E.  Parker,  and  C.K.  Birdsall,  "Particle  Transport  due  to  Kelvin-Helmholtz  Vortices  and  Small 
Scale  Turbulence." 

US-Japan  Workshop  on  "Advanced  Computer  Simulation  Techniques  Applied  to  Plasmas  and  Fusion," 
UCLA,  September  26-28,  1990. 

V.  Vahedi,  J.P.  Verboncoeur,  M.  Surendra,  and  C.K.  Birdsall,  "A  Monte-Carlo  Collision  Model 
for  the  Particle-in-Cell  Method." 

Course 


M.A.  Lieberman  and  C.K.  Birdsall,  "Low  Pressure  Plasma  Discharges  Used  in  Semiconductor  Process¬ 
ing",  UC  Extension,  San  Francisco,  CA,  May  23-25  1990,  and  Oxford,  England,  July  9-11,  1990. 


Invited  Talks 


C.K.  Birdsall,  "Interactive  Plasma  Computer  Experiments;  Plasma  Device  Simulations  on  PC’s 
and  Workstations,"  January  16,  1990  at  Naval  Research  Laboratory,  Washington,  D.C. 

C.K.  Birdsall,  "Magnetized  plasma  sheath,  with  Kelvin-Helmholtz  instability,  vortices,  Bohm-like 
diffusion",  and  "Computer  experiments  with  bounded  collisional  plasmas  on  fast  PC’s,  March  26- 
27,  1990,  at  University  of  Wisconsin,  Madison,  Wisconsin. 

C.K.  Birdsall,  "Interactive  Plasma  Computer  Experiments;  Plasma  Device  Simulations  on  PC’s 
and  Workstations,"  June  1, 1990  Naval  Post  Graduate  School,  Monterey,  CA. 

C.K.  Birdsall,  "Particle-In-Cell  Combined  with  Monte  Carlo  Collisions  In  Living  Color:,  October 
16-19, 1990,  at  the  43rd  Annual  Gaseous  Electronics  Conference,  Urbana-Champaign,  Illinois. 


Presented  at  Sherwood  Fusion  Theory  Conference 
Williamsburg,  VA,  April  23-25,  1990 


Particle  Simulations  Of  Transport  In  A  Diverted  Tokamak 
Scrape-Off  Layer:  The  High-Recycling  Regime* 

Richard  J.  Procassini  and  Charles  1C  Birdsall 

Electronics  Research  Laboratory 
University  of  California 
Berkeley,  CA  94720 

Bruce  I.  Cohen 

Magnetic  Fusion  Energy  Division 
Lawrence  Livermore  National  Laboratory 
Livermore,  CA  94550 

The  transport  of  partides  and  energy  through  a  tokamak  scrape-off  layer  (SOL)  is  stucfied  via  the  particle 
simulator)  code  DIPSI  (Direct  Implicit  Plasma  Surface  Interactions).  This  code  combines  direct-implicit 
particle-in-cell  (PIC)  techniques  with  Monte  Carlo  models  for  Coulomb  collisions,  and  charged/neutral  atomic 
physics  interactions.  The  fully-kinetic  guiding-center  PIC  code  provides  a  self-consistent  solution  of  the 
electrostatic  potential  profile  in  one  spatial  dimension  (along  the  open  field  lines  in  the  SOL)  and  two  velocity 
components  (vvJ< 

The  Monte  Carlo,  binary-partWe  Coulomb  colston  model  conserves  both  the  momentum  and  kinetic  energy 
of  the  interacting  particles.  The  fufl  range  of  Coulomb  coBsionai  processes  (e-e,  e-i  and  H)  are  Included. 
Recycled  neutral  particles  are  treated  in  an  ad  hoc  manner,  the  density  is  assumed  to  decay  exponentially 
with  increasing  dstance  from  the  dtartor  plate  (the  e-fokfing  length  is  taken  to  be  the  average  mean-free 
path  to  iontzatton)  and  the  temperature  of  the  neutrals  is  prescribed  by  the  user.  The  neutral  particles  interact 
with  the  background  plasma  via  charge-exchange  events  with  ions  and  impact  ionization  events  arising  from 
electron/neutral  coiGsions. 

This  model  is  used  to  determine  the  effect  of  charged/neutral  interactions  on  the  transport  properties  of  a 
tokamak  SOL  plasma.  The  prescribed  neutral  particle  density  wBI  span  a  wide  range  of  recycling  fractions. 
Values  of  the  presheath  and  collector  sheath  potential  drops,  plasma  temperature,  flow  velocity  and  parallel 
heat  fluxes  as  a  function  of  charged/neutral  cottsionality  wll  be  presented. 

We  wfll  also  dscuss  the  possfeiHty  of  extencSng  this  model  to  include  the  effects  of  sputtered  impurity  particles, 
which  could  be  evolved  as  a  second  ionic  species.  This  would  allow  one  to  determine  the  effect  of  sputtered 
impurities  on  the  potential  profBe  (and  uRimately ,  the  transport  of  partides  and  energy)  in  the  plasma-sheath 
region. 


*  Work  performed  for  the  U.S.  Department  of  Energy  by  the  Lawrence  Livermore  National  Laboratory  under 
contract  W-7405- ENG-48. 


Presented  at  IEEE  International  Conference  on 
Plasma  Science 
Oakland ,  CA 
May  21-23,  1990 
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Particle  Simulations  of  Transport  in  a  . 
High-Recycling  Divertor  Scrape-Off  Layer 

Richard  J.  Procassini  and  Charles  K.  Birdsatl 

Electronics  Research  Laboratory 
University  ot  California 
Berkeley,  CA  94720  U.S.A 

Bruce  I.  Cohen 

Magnetic  Fusion  Energy  Division 
Lawrence  Livermore  National  Laboratory 
Livermore.  CA  94550  U.S.A. 

The  transport  ot  particles  and  energy  through  a  tokamak  scrape-off  layer 
(SOL)  is  studied  via  the  particle  simulaton  code  DIPS!  (Direct  Implic* 
Plasma  Surface  Interactions).  This  code  combines  direct-implicit 
particle- in -cell  (PIC)  techniques  with  Monte  Carlo  models  lor  Coulomb 
collisions  and  charged/neutral  atomic  physics  interactions .  The  fully -kinetic 
guiding-center  PIC  code  provides  a  sell -consistent  solution  of  the  elec¬ 
trostatic  potential  profile  in  one  spatial  dimension  (along  the  open  field  Knee 
in  the  SOL)  and  two  velocity  components  (*„v,) . 

The  Monte  Carlo,  binary-particle  Coulomb  collision  model  conserves  both 
the  momentum  and  kinetic  energy  of  the  interacting  particles.  The  lu#  range 
of  Coulomb  coUisional  processes  (e-e.  e-i  and  M)  are  included.  Recycled 
neutral  particles  are  treated  In  an  ad  hoc  manner  the  density  is  assumed 
to  decay  exponential  with  increasing  distance  from  the  divertor  plate  (the 
e-folding  length  is  taken  to  be  the  average  mean-free  path  to  ionization) 
and  the  temperature  of  the  neutrals  is  prescribed  by  the  user.  The  neutral 
particles  interact  with  the  background  plasma  via  charge-exchange  events 
with  ions  and  impact  Ionization  events  arising  from  electron/neutral  ces¬ 
sions. 

This  model  is  used  to  determine  the  effect  ot  chargad/heutral  interactions 
on  the  transport  properties  of  a  tokamak  SOL  plasma.  The  prescribed 
neut  ral  particle  density  wfll  span  a  wide  range  of  recycling  fractions.  Values 
of  the  presheath  and  cofcctor  sheath  potential  drops,  plasma  temperature, 
flow  velocity  and  parallel  heat  fluxes  as  a  function  of  charged/neutrai  col- 
HsionaWy  will  be  presented. 

Extensions  of  this  model  which  include  the  effects  of  sputtered  impurity 
particles  w*  also  be  dtocussed.  The  sputtered  particles  may  be  evolved 
as  a  second  ionic  spedee.  This  wM  alow  one  to  determine  the  effect  of 
sputtered  impurities  on  the  potential  profile  (and  ultimately,  the  transport  of 
particles  and  energy)  in  the  plasma-sheath  region. 
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SHEATH  VOLTAGE  RATIO  FOR 
ASYMMETRIC  RF  DISCHARGES' 

M.  V.  Alves2,  M.  A.  Lieberman, 

V.  Vahedi  and  C.  K.  Birdsall 
University  of  California 
Berkeley,  CA  94720 

Spherical  and  cylindrical  many-particle  models  are  being 
used  to  simulate  RF  discharges  in  which  the  RF  powered  and 
the  grounded  electrodes  have  different  areas.  This  asymmetry 
determines  the  magnitude  of  the  self-bias  voltage  V ,  (the  ion 
bombarding  energy)  at  the  powered  electrode,  which  is  a 
critical  process  parameter.  Recent  analytical  models  for  the 
sheath  voltage  ratio  including  the  effect  of  the  floating 
potential  have  been  developed1  4  that  agree  with  some 
experimental  results.  Onedimensional(radial)  spherical  shell 
models  have  also  been  developed1,  incorporating  various 
assumptions  for  the  sheath  and  the  glow  discharges,  leading 
to  a  scaling  which  is  in  agreement  with  some  measurements. 
We  have  simulated  the  spherical  model  with  a  non  uniform 
ionization  and  the  results  agree  w  ith  the  theory.  The  cylin¬ 
drical  simulation  shows  that  the  floating  potential  plays  an 
important  role.  The  simulation  results  will  be  given 
graphically  as  the  sheath  voltage  ratio  versus  area  ratio,  and 
will  be  compared  with  the  theory.  The  simulation  codes  are 
PDC1  (cylindrical)  and  PDS1  (spherical)  which  utilize 
particle-in-cell  techniques  plus  Monte-Carlo  simulation  of 
electron-neutral  (elastic,  excitation  and  ionization)  and  ion- 
neutral  (scattering  and  charge-exchange)  collisions6  7. 


1.  This  work  is  supported  in  pan  by  ONR  contract  N00014- 
90-J-1198. 

2.  M.  V.  Alves,  visiting  from  INPE-Brazil,  supported  par¬ 
tially  by  CAPES-Ministry  of  Education,  Brazil. 

3.  M.  A.  Lieberman,  to  appear  in  J.  Vac.  Sci.  Tcchnol.  A , 
(1990) 

4.  A.  M.  Pointu,  Appl.  Phys.  Lett.,  50:1047,  (1987) 

5.  M.  A.  Lieberman,  J.  Appl.  Phys.,  65:4186,  (1989) 

6.  M.  V.  Alves,  V.  Vahedi  and  C.  K.  Birdsall,  Bull.  APS, 
34:2028,  (1989)  (Abstract) 

7.  Codes  available  from  Industrial  Liaison  Program,  EECS 
Dept.,  UC  Berkeley. 
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WORK  DONE  AND  ENERGY  BALANCE 
IN  RF  DISCHARGES1 

J.  P.  Verboncoeur  and  V.  Vahedi 
M.  A.  Liebetinan,  and  C.  K.  Biidsall 
University  of  California 
Berkeley,  CA  94720 

The  fields  applied  to  a  parallel  plate  RF  discharge  do  work 
on  the  plasma  panicles  (deposit  energy)  non-uniformly  in 
space  and  time  at  a  rate  given  by  J  •  E.  Similarly,  energy  is 
lost  by  the  plasma  particles  when  they  strike  the  electrodes 
or  collide  with  neutrals  or  radiate.  We  are  currently  charac¬ 
terizing  the  work  done  and  the  losses  as  a  function  of  time 
and  space  in  a  manner  similar  to  that  of  Vender  and  Boswell2 3 
self-consistently.  We  compare  these  results  to  discharge 
models  and  other  simulations.  The  simulation  code  is  PDP1 
which  utilizes  panicle-in-cell  techniques  plus  Monte-Carlo 
simulation  of  electron-neutral  (elastic,  excitation  and  ion¬ 
ization)  and  ion-neutral  (scattering  and  charge-exchange) 
collisions’ 4.  The  code  can  simulate  different  environments 
with  various  gases  and  various  pressures  commonly  used  in 
RF  discharges  needed  in  plasma  processing. 


1.  This  work  is  supported  in  part  by  ONR  contract 
N00014-90-J-1198 

2.  D.  Vender  and  R.  W.  Boswell.  Measurements  of  Power 
Dissipation  in  a  RF  Plasma  Simulation.  Report  ANU-PRL- 
PP88/8,  Plasma  Research  Laboratory,  Australian  National 
University,  Canberra  (1988). 

3.  1.  J.  Morey,  V.  Vahedi  and  J.  Verboncoeur,  Particle 
Simulation  Code  for  Modeling  Processing  Plasmas,  fiulL 
AES.  34:2028.  (1989)  (Abstract) 

4.  Codes  available  from  Industrial  Liaison  Program.  EECS 
Dept.,  UC  Berkeley. 
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A  COLLISIONAL  MODEL  FOR  PLASMA 
IMMERSION  ION  IMPLANTATION1 

V.  Vahedi,  M.  A.  Lieberman, 

M.  V.  Alves2,  J.  P.  Verboncoeur, 
and  C.  K.  Biidsall 
University  of  California 
Berkeley,  CA  94720 

In  plasma  immersion  ion  implantation,  a  target  is  immersed 
in  a  plasma  and  a  series  of  negative  short  pulses  are  applied 
to  it  to  implant  the  ions.  A  new  analytical  model  is  being 
developed  for  the  high  pressure  regimes  in  which  the  motion 
of  the  ions  is  highly  collisions!  The  model  provides  values 
for  ion  flux,  average  ion  velocity  at  the  target,  and  sheath  edge 
motion  as  a  function  of  time.  The  model  suggests  that  the 
transient  ion  flux  at  the  target  scales  with: 

> 

Jc  aVfeuortoQ+av)  ’ 

where  tUg  and  14,  are  characteristic  frequency  and  velocity  of 
the  ions  in  the  sheath,  rig  is  the  density,  and  r  is  the  time.  These 
values  are  being  compared  with  those  obtained  from  simu¬ 
lation  and  show  good  agreement  A  review  will  also  be  given 
(for  comparison)  of  the  earlier  work  done  at  low  pressures, 
where  the  motion  of  ions  in  the  sheath  is  collisionless’,  also 
showing  good  agreement  between  analysis  and  simulation4. 
The  simulation  code  is  PDP1  which  utilizes  particle-in-cell 
techniques  plus  Monte-Carlo  simulation  of  electron-neutral 
(elastic,  excitation  and  ionization)  and  ion-neutral  (scattering 
and  charge-exchange)  collisions’. 


1.  This  work  is  supported  in  pan  by  ONR  contract 
N00014-90-M198 

2.  M.  V.  Alves,  visiting  from  INPE-Brazil,  supponed  par¬ 
tially  by  CAPES-Ministry  of  Education,  Brazil. 

3.  M.  A.  Lieberman,  J.  Appl.  Phys.,  66:2926,  (1989) 


4.  I.  J.  Morey,  V.  Vahedi  and  J.  Verboncoeur,  Particle 
Simulation  Code  for  Modeling  Processing  Plasmas,  Bui! 
AES,  34:2028.  (1989)  (Abstract) 

5.  Codes  available  from  Industrial  Liaison  Program,  EECS 
Dept.,  UC  Berkeley. 
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Particle  Transport  due  to  Kelvin-Helmholtz  Vor¬ 
tices  and  Small  Scale  Turbulence  :  S.  E.  Parker  and 
C.  K.  BlRDSALL,  University  of  California,  Berkeley — In  simulations 
of  the  cross-field  plasma  sheath,  which  has  B  parallel  to  an  absorb¬ 
ing  wall  and  U||  =  fcj|  =  0,  we  have  found  that  the  E  x  B  drift  velocity 
shear  produces  the  Kelvin-Helmholtz  instability  which  saturates  to  large 
scale  drifting  vortices  as  -2,  diameter  as  5-l5pt),  with  subsequent 
Bohm-like  diffusion.1  In  addition  to  the  circular  vortex  flow,  there  is 
a  smaller  amplitude,  small  scale  turbulent  spectrum  ( <  1,  kpi  ~  1 
and  u>  —  uici).  We  have  now  modeled  the  vortex  as  a  Gaussian  potential 
well  and  the  small  scale  turbulence  as  perturbing  plane  waves  traveling 
parallel  to  the  wall.  With  no  waves  present,  E  X  B  drifting  test  elec¬ 
trons  stay  on  the  equipotential  lines.  The  addition  of  one  wave  leads  to 
stochastic  motion  and  loss  of  particles.  Use  of  many  waves  of  a  given 
spectrum  will  also  be  presented,  as  well  as  ties  between  rate  of  trans¬ 
port  and  wave(s)  amplitude.  Our  objective  is  to  understand  how  the 
K-H  vortices  provide  a  mechanism  for  enhanced  transport. 

'This  work  supported  by  DoE  contract  No.  DE-FG03-S6ER53220  and 
ONR  contract  N00014-85-K-0809. 

XK.  Theilhaber  and  C.  K.  Birdsall,  Phys.  Fluids  B  1,  2241,  Parts  I  and 
II,  (1989). 
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